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Abstract. Recent theoretical developments concerning Higgs production at the Large Hadron Collider 
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1. Introduction 



The upcoming searches for and anticipated studies of Higgs bosons at the Large Hadron CoUider 
(LHC) have been the driving force for interesting theoretical developments for many years. Quite 
a number of observables and input quantities to experimental analyses are meanwhile known to 
high precision. In particular, the most important Higgs production cross sections are under good 
theoretical control: the next-to-leading order (NLO) predictions for associated tiH production [1-3] 
and weak boson fusion [4], as well as the next-to- next-to-leading order (NNLO) cross sections for 
gluon fusion [5-7] and Higgs Strahlung [8,9] all exhibit nicely converging perturbative series and 
only a moderate dependence on the renormalization and factorization scale. 

This brief overview should be considered as a report on recent developments concerning higher 
order calculations for neutral Higgs production at the LHC. For more comprehensive surveys, we 
refer the reader to some recent reviews (see, e.g., Refs. [10,11]). 



2. Gluon fusion 



Higher order corrections to the gluon fusion process have been of great interest for many years now. 
The main reason is, of course, that it is one of the most important discovery channels for Higgs 
bosons at the LHC, and that therefore its cross section has to be under good theoretical control. In 
fact, the NLO radiative corrections turned out to be very large [12-14], amounting to an increase 
of up to 100% with respect to the leading-order (LO) value. In addition, the NLO corrections did 
not lead to a decrease of the renormalization and factorization scale dependence (when measured 
in absolute rather than relative values of the cross section). Due to these issues in the theoretical 
prediction, the if-factor was often neglected in experimental analyses. 

A lesson learned from the NLO corrections was that the gluon-Higgs interaction seems to be 
approximated very well [15,16] by an effective Lagrangian 



1 



if the LO top mass dependence of the cross section is factored out. In Eq. (1), v = 246 GeV 

and C{as) is the Wilson coefficient which is meanwhile known through [17,18]. This obser- 
vation allowed to tackle the NNLO calculation on the basis of Eq. (1) [19-21]. The fuU NNLO 
corrections [5-7] exhibited the features of a well-behaved pcrturbative series: they are significantly 
smaller than the NLO corrections, and also the scale dependence reduces to an acceptable level. 

Progress concerning higher order corrections to the gluon fusion process has been made in various 
respects. On the one hand, the validity of the perturbative prediction for the total cross section 
has been confirmed by the evaluation of effects that go beyond NNLO. On the; other hand, various 
resummations for kinematical distributions of the Higgs boson have been carried out. And finally, 
a fully differential partonic Monte Carlo program, valid through NNLO, has been developed. Let 
us discuss these topics in more detail in what follows. 

Inclusive Higgs production. With the higher order QCD effects for gluon fusion being quite sizable, 
one may wonder about the reliability of the fixed order NNLO prediction. To answer this question, 
one may try to identify and resum the dominant terms to the inclusive cross section. In fact, it had 
been realized long ago that soft gluon radiation contributes significantly to the total rate [15,20,21]. 
However, how well the cross section is approximated by this contribution alone depends strongly 
on the way the "soft limit" is defined. To see what wc mean by this, consider the general expression 
for the hadronic cross section a{s) in terms of parton densities {i = q, q, g) and the partonic 
cross section a{s): 

cr(s) = / da;i / diX2<i)i{xi)4>j{x2)aij{xiX2s) . (2) 
Jo Jo 

One way to define the soft limit is to expand a in the limit x = Mfj/s — > 1: 

aij{s)=ao\a5{l-x)+^bkVk{x) + --- I (3) 



where 

Vk{x) 



k>0 



ln''(l-x) 
1-x 
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and the dots denote formally subleading terms which are dropped. However, one may equally well 
trade a factor of x between the partonic cross section and the parton density functions, and rather 
expand a{s)/x around x = 1. This will lead to significantly different numerical results for the 

hadronic cross section a [20,21]. 

A definition of the soft limit that seems to approximate the full NNLO cross section very well has 
been found in Ref. [22]: one transforms a{s)/x from x space to Mellin moment space ('W space"), 
and drops all terms that are of order 1/N. The terms a^ln*^ N were then resummed to next-to- 
next-to-leading logarithmic accuracy (NNLL), with a rather moderate numerical impact [22]. This 
indeed indicates the aforementioned stability of the NNLO expression. 

The most recent achievement concerning the fixed-order calculation is the evaluation of the soft 
terms through N^LO or, in terms of Eq. (3), the coefficients bk, fc = 0, . . . , 5, through [23]. The 
(5(1 — a;) piece in Eq. (3) receives contributions from the virtual terms which are still unknown, but 
one may deduce from the lower order results that they are numerically small. Note also that the 
terms for k = 1, . . . ,5 could be derived from the NNLL resummation formula [22], thus providing 
a useful check. The effect of these N^LO terms is again a mild increase of the cross section (for 
I^F = I^R ~ Mh), and a reduction of the scale uncertainty. On the other hand, these terms allow to 
push the resummation of the soft terms to higher orders [23,24], with again rather small numerical 
impact (see Fig. 1). 




It remains to be said that apart from the phenomenological significance of such higher order results, 
they appear to reveal interesting structures of the perturbative series that are worth studying for 
their own sake [25,26]. 

Clearly, if so much effort is put into minimizing the theoretical uncertainty due to QCD effects, 
one needs to start thinking about electro-weak corrections as well. In fact, terms of order Gprn^ 
have been known for quite a while now, resulting in effects below 1% of the LO rate [27]. The 
full set of Feynman diagrams can be divided into those that contain a top quark, and those that 
contain only light quarks. The latter set has been evaluated in Ref. [28], while the full result for 
the former was obtained in Ref. [29] . The numerical impact of the electro- weak corrections ranges 
between 5 and 8% of the LO term. 

Transverse momentum and rapidity distributions. A large ii'-factor for the total cross section 
leaves open the question on how the radiative corrections affect different regions of phase space. 
Better insight into this issue is provided by differential quantities such as dcr/dpTdy, where pr and 
y are the transverse momentum and the rapidity of the Higgs boson, respectively. 

At the partonic level, the Higgs boson can only be produced at finite transverse momentum px 
if the latter is balanced by the real radiation of a quark or a gluon. Thus, distributions at non- 
zero pt are related to the process i?+jet whose LO prediction is of order a^. The corresponding 
NLO effects were studied both analytically [30,31], and in the form of a partonic Monte Carlo 
program [32], yielding a rather flat dependence of the /^-factor on pT and y at intermediate values 
of these variables. At small p^, the fixed-order perturbative approach breaks down, but this can be 
accounted for by resummation of logarithms (see, e.g., Ref. [33] and references therein). At large 
Pt, on the other hand, one again encounters similar logarithms as for the inclusive rate, arising 
from soft gluon radiation. Their resummation is known through NLL [34] . 

NNLO Monte Carlo. Currently the most general higher order prediction for the gluon fusion 
process is available in the form of the partonic NNLO Monte Carlo program FEHiP [35]. It allows 
to study arbitrary kinematical distributions of the Higgs boson with NNLO accuracy, as well as the 
application of phase space cuts. This provides detailed information on how the radiative corrections 
affect the various regions of phase space. 

For example. Table 1 compares the ratio A'^^^ of the NNLO to the NLO cross section as obtained 
by the fully inclusive calculation (subscript "inc") to the one where "standard cuts" are applied 
(subscript "cut"; see Ref. [35] for details). The numbers show that the radiative corrections are 
only slightly affected by the cuts, and that the cross section including cuts is approximated to 
better than 5% by the quantity K^"^^^ x o'n'lo- 
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Table 1. Comparisons between the cut and inclusive cross sections for dif- 
ferent Higgs masses. The second column contains the ratio of the NNLO 
cross section with the standard cuts over the inclusive cross section, while the 
third column contains the ratio of cut and inclusive results for the /f-factor 
K^'^^ = aNNLo/o-NLO- It is /iR = /iF = Mr 12. From Ref. [35]. 
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This observation motivates another step towards more reaUstic higher order event simulations: In 

order to transfer the purely partonic NNLO result of Rcf. [35] to truly hadronic final states, the 
Pythia [36] and MCONLO [37,38] event generators have been supplemented by a re-weighting grid 
in the pr-y plane [39], evaluated from the ratio of the partonic FEHiP and the hadronic results, 
integrated over two-dimensional intervals. This procedure is based on the fact that sufficiently 
inclusive quantities should be described equally well in a partonic and a hadronic approach. A 
similar strategy was followed in Ref. [40], where the re- weighting was based only on the NLO 
PT-spectrum of the Higgs boson, however. 

Background calculations. An important issue for Higgs searches and studies at the LHC is the 
theoretical control of background processes. The number of higher order results available in this 

context is way too large to even attempt giving proper credit to each one of them. An extensive 
list of programs to evaluate higher order cross sections can be found at Rcf. [41]. 

Quite often, side band subtractions rather than NLO simulations will be the most efficient way in 
order to separate the signal from the background, once enough data are available. But in certain 
cases, such a procedure will not be possible, for instance if missing energy in the Higgs decay 
does not allow the reconstruction of a Higgs mass peak. An example for such a case is the WW 
decay mode of the Higgs boson. In fact, in order to use this channel as a discovery mode, one 
needs to keep track of the angular correlations among the Higgs decay products [42] . Also here, 
the NLO corrections have been known for a while. However, recently it was found that the gg 
initiated component, although being formally of NNLO, can amount to 30% of the NLO rate, once 
the relevant cuts for Higgs searches are applied [43,44]. 



3. Weak boson fusion 

The weak boson fusion (WBF) process itself is under very good theoretical control: the NLO cor- 
rections have a comparatively simple structure, since single gluon exchange between the incoming 
quarks is not allowed by color conservation, meaning that 5-point functions are absent in the cal- 
culation. The NLO corrections are available in the form of a partonic Monte Carlo program [4], 
allowing for application of cuts as it is particulary important for this process in order to separate 
it from the background. 

The challenge concerning radiative corrections is indeed related to these background processes. 
The dominant source is Higgs production in gluon fusion, when the Higgs is associated with two 
jets. The LO prediction is of order , meaning that the renormalization scale dependence is rather 
large. The full top mass dependence of the cross section at LO has been evaluated in Ref. [45]. 
The NLO calculation involves massive two-loop five-point functions and is certainly out of reach 
at the moment. However, the LO calculation revealed that for jet transverse momenta pTj ^ mi, 
one may integrate out the top quark, thus arriving at one-loop five-point functions, calculated in 
Ref. [46]. However, these still need to be supplemented by the real radiation contribution, the 
amplitude for which has been evaluated in Ref. [47-49] . 

Other important background processes to WBF are Vjj and VVjj production, and also here NLO 
corrections are available [50,51]. 



4. Supersymmetry 

Gluon fusion remains one of the most important production modes also in supersymmetric models. 
In fact, since the pseudo-scalar Higgs boson in the MSSM has no tree-level coupling to vector 
bosons, it cannot be produced through WBF, for example. Therefore, gluon fusion and associated 
bbA production are the dominant production modes in this case (see, e.g., Ref. [52]). 

Many of the higher order results that are available for SM Higgs production can be taken over 
to the case of the neutral, CP-even Higgs bosons within the MSSM. For example, if the squarks 



are heavy, they do not contribute significantly to the gluon-Higgs couphng, which is then again 
mediated predominantly by top and, for not too small values of tan/?, bottom loops. The cross 
section for h, i/-production including QCD corrections can then be derived easily from the SM 
expression. 

Pseudo-scalar Higgs production, on the other hand, requires a modification of the top-Higgs and 
thus the effective gluon-Higgs coupling with respect to Eq. (1). But also here, many higher order 
corrections are known, for example the inclusive NNLO cross section [53,54,7], and various NLO 
distributions [55-57]. 

Larger values of tan/3 increase the importance of bottom loops to the gluon-Higgs coupling (see, 
e.g., Refs. [16,58]). An effective theory for the gluon-Higgs interaction along the lines of Eq. (1) is 
not known in this case, such that higher order calculations are much more difficult. In fact, only the 
NLO result is available at the moment, in the form of a one-dimensional integral representation [59]. 
Let us remark, however, that the virtual corrections have meanwhile been expressed in terms of 
analytic functions [60]. 

If their masses are not too large, top squarks may influence the gluon-Higgs coupling as well. In 
this case, the NLO gets more involved as compared to the SM, mostly because several mass scales 
enter the problem. However, one may again employ an effective theory approach analogous to 
Eq. (1), where top quarks and squarks are considered as heavy [61,62]. The particle spectrum of 
the effective theory is then the same as in the SM case, and the only unknown quantity is the Wilson 
coefficient C(as). The latter is obtained from massive tadpole integrals which can be evaluated 
analytically through two loops using the proper reduction formulas [63] . 

This allows one to evaluate the NLO corrections to SUSY Higgs production for small values of tan /3 
(for large tan/3, the bottom and sbottom loop effects may be taken into account approximately 
using the leading order expression and resummation of tan/3 terms [64]). The results for both the 
production of a CP-even [62] and a CP-odd [65] Higgs boson have been evaluated through NLO 
in this way. For the CP-even Higgs, even an estimate at NNLO has been obtained [66] (denoted 
NNLO' in Fig. 2 below). 

A particulary dramatic scenario is given by the so-called "gluophobic Higgs" [67,68], where the 
quark and squark loops interfere distructively, such that the Higgs coupling to gluons becomes very 
small. Fig. 2 shows the effects of higher orders in as in this region of SUSY parameter space. One 
observes that the radiative corrections do not change the general behaviour of the cross section. 
Rather, the LO rate is multiplied by an almost constant i^-factor close to the one of the SM 
calculation [62]. 
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Figure 2. Thick lines: gluon fusion cross section in the MSSM (for details, 
see Ref. [62]). The thin lines show the result when the squark effects are 
neglected. 




Concerning differential distributions, the LO diagrams to Higgs plus jet production due to squark 
loops have been evaluated in Refs. [69,57]. For higher order effects, one currently needs to rely on 
the effective Lagrangian approach due to the complexity of the calculation. Then, however, the 
SUSY effects factorize into the Wilson coefficient, just like for the inclusive rate. 



5. Bottom quark annihilation 



In SUSY, Higgs production in association with bottom quarks can give a significant contribution 
to the total Higgs production cross section at the LHC. In fact, it can even exceed the gluon fusion 
component. The proper theoretical description has been a subject of discussion for quite some 
time now. The difficulties arise from the fact that a potentially large mass difference between the 
Higgs boson and the bottom quark can lead to large logarithms that originate from integration 
over the collinear region of one or both of the produced bottom quarks. It has been suggested to 
use bottom quark parton densities as a way to resum these logarithms. However, this so-called 
5-flavor scheme (5-FS) and the 4-flavor scheme (4-FS), where these logarithms are not resummed, 
lead to considerably different numerical results for the total cross section. 

It was then realized that the discrepancies between the two approaches are much smaller once the 
factorization scale /ip for the bottom densities is chosen significantly lower than the supposedly 
"natural" choice fip = Mh- In fact, based on the argument that factorization works only in the 
collinear limit, it was suggested that a reasonable choice was //p — Mh/^ in this case [70-72]. This 
is because for px > the pT-distribution of the bottom quarks in the final state begins to 

deviate significantly from the collinear form Aa/dpT ^ 1/pt [73]. 

This choice for the factorization scale later has found support from the NNLO result for that 
process, evaluated in the 5-FS [74]. This is shown in Fig. 3 (a) for the LHC. Clearly, convergence 
of the perturbative series appears to be much better for scales below Mh rather than above. 

Yet another intriguing feature of this NNLO result indicates that indeed the scale ^p = M^f/4 is 
markedly different from any other. To see this. Fig. 3 (b) shows separately the contributions from 
the bb, bg + bg, gg, and other, much smaller partonic sub-processes in the MS scheme. Also shown 
is the sum of all these curves (solid line), corresponding to the NNLO result. It so happens that 
right at /ip = M^f/4, all contributions except for the bb term practically vanish simultaneously [75]. 

One may wonder whether the bottom quarks can be taken massless for the gg component of the 
NNLO contribution, which does not have an initial state bottom quark. These effects, however, 
are expected to be of order m^/Mfj and thus negligible. This is indeed observed when comparing 
the massive [71] to the massless [74] result of this component [75]. 
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Figure 3. (a) LO, NLO, and NNLO prediction for the inclusive bbH cross 
section in the 5-FS [74] . (b) Individual components to the NNLO prediction 
arising from various subprocesses in the MS scheme. See also Ref. [75] 
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Clearly, the pure 5-FS can not be applied directly to exclusive bbH production, i.e., if one or both 
of the bottom jets arc required to be produced at large transverse momenta. In the fully exclusive 
case, the NLO corrections are available [76,77]. If only one bottom quark is required at large pt, 
one may apply the 5-FS for the one that is integrated over. Also in this case, the NLO corrections 
are known [78,79]. 

Concluding this section, the bbH process has been a very inspiring subject over the past few years, 
and there would still be several aspects to be discussed. However, in this short write-up, we have 
to refer the interested reader to the recent literature (see, e.g., Refs. [75,80,81] and references 
therein) . 



6. Conclusions 

Higgs physics has been a very fruitful field of research and, with the LHC data to come, will be 
even more so in the future. Progress has been fast-paced, and many of the results and techniques 
are general enough to find applications also in very different contexts. I have tried to summarize 
the most significant developments of the past few years related to Higgs production at the LHC, 
and to direct the reader to the relevant literature whenever more detailed information is required. 
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